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Abstract 

Dit verslag belichaamt het onderzoek dat werd verwezenlijkt in kader van ons 
bachelor project. Het onderzoek bestond eruit om na te gaan of er een cor¬ 
relatie bestaat tussen de orbitale periode van het binair systeem van Barium 
sterren en hun barium abundantie. Deze correlatie werd onderzocht aan de 
hand van een verzameling van 9 Barium sterren, hierbij werden de belangrijk¬ 
ste stellaire parameters geverifieerd door middel van het onderzoeken van de 
stellaire spectra. Vervolgens werd de barium abundantie bepaald aan de hand 
van de absorptielijnen corresponderend met het element barium. Op basis van 
deze resultaten kon geconcludeerd worden dat er geen significante correlatie 
bestaat tussen beide parameters. 


1 Introduction 

In 1835 Auguste Comte, a French philosopher, predicted that the Chemical com- 
position of the Sun or any other star would never be known for the simple rea- 
son that mankind cannot reach them. Some 30 years later, Bunsen and Kirch- 
hoff had made a first step to prove Comte was wrong by observing the spectra 
of hot gases and by comparing them with the spectrum of the Sun. It became 
clear that the Sun and any other star was made out of a mix of hot gasses 
sending out electromagnetic radiation. They found that by analyzing the elec- 
tromagnetic radiation, it is possible to obtain the main characteristics about an 
observed star. 


1.1 Barium Stars 


In 1951 William P. Bidelman and Philip C. Keenan discovered a new group of 
special G- and K-type red giants, characterized by strong spectral lines of ion- 
ized barium and certain carbon rich molecules like CH, C 2 and CN |(Bidelman| 
|& Keenan, 1951) The fact that those BA II -lines (first excitation of barium) 
were so strong was abnormal and therefore this group of stars was labelled as 
"Barium stars". They became a hot subject in the scientific world because of 
their peculiar Chemical composition, even more when Barium stars were dis- 
covered on the main sequence of the Hertzsprung-Russell (H-R) diagram. At 
this evolutionary stage the abnormal abundance of barium cannot be explained 
by nucleosynthesis nor as an element present during the formation of the star. 
Hence, this peculiar abundance of Barium can only be caused by external fac¬ 
tors. 


Robert D. McClure showed in 1980 that each Barium star is a member of a 
binary system accompanied by a white dwarf (McClure et al., 1980)| He de- 
fined the Barium stars as Population I red giants of spectral type G-K having 
enhanced abundances of carbon and s-process elements. The s-process or the 
slow neutron capture is the process responsible for the nucleosynthesis in AGB 
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(Asymptotic giant branch) stars. 

McClure also had an hypothesis to explain the enhanced barium abun- 
dances. In a binary star one star will be more massive than the other. The more 
massive star will evolve faster than its companion. When the more massive 
star reaches the AGB, its radius is much larger than initially This can induce 
a mass transfer because the outer layers of the star are far away from the star 
itself and starts to feel the gravitational attraction of the companion. After the 
mass transfer the star becomes a white dwarf and stops evolving. But the other 
star that hasn't reach the AGB-stage yet has already a high barium abundance. 
Therefore mass transfer is the main reason why the less massive star gets pol- 
luted with barium and other s-process elements. On the figure below ([TJ the 
evolution is schematically shown for a binary system where one of the stars 
pollutes the other. 
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Figure 1: Evolution of a binary system of two stars to a system with a Barium 
star and a white dwarf 


KJorissen et al., 1996) 


1.2 Binary Systems 

The stars we observe at the night sky can be divided into two categories. There 
are the truly single stars, as for example our Sun, and there are pairs of stars, 
with two (or more) components moving in bound orbits about their barycen- 
tre, which we call binary systems of stars, or simply binary stars. Since this 
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bachelor project is concerned about Barium stars it is clear that our may focus 
will be on this second category of stars, namely binary stars. Now, there are 
several types of binary systems. 

The first binary systems discovered were the visual binaries in which case, 
their relative motions can be resolved into two components, and their apparent 
motions on the night sky can be measured. Another very distinguishable type 
of binary systems is the eclipsing binary When the plane in which the two 
stars rotate around their barycentre is almost parallel with the line of sight 
of the observer, a difference in luminosity can be noticed when both stars are 
eclipsing. As a result, this difference in luminosity allows the determination of 
the radial velocity of the star. 

A third type of binary system is the astrometric binary star. At first sight 
some stars have weird apparent motions on the sky. This unnatural path can 
for certain stars be explained by the presence of a companion. The presence 
of this companion affects the path of the main bright star due to gravitational 
interaction. 



Tangenüal motton 
ns dopplof shffi 


ïangeriMl motton 
no ttoppler shtH 


A Spectroscopie Binary System 
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Figure 2: Example of a spectroscopie binary system 
[Spectroscopie Binariës] 
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The last type of binary stars is the spectroscopie binary As a result of the 
relative motion of the star with respect to the observer a Doppler shift can occur 
in the spectrum of the star. If the relative motion of the star is towards the 
observer, its spectrum is shifted to the blue. Inversely if the relative motion is 
away from the observer, the spectrum is shifted to the red. In a binary system 
the gravitational interaction causes the stars to rotate around the barycentre, 
this means that the spectra of the stars in a binary system will be periodically 
shifted to the red and to the blue. This means that if a spectrum is alternatively 
shifted to the red and to the blue, the star might be a part of a binary system. 
An example of a spectroscopie binary is given in Figure [2] The binaries that 
will be discussed in this project are all spectroscopie binaries. 

1.3 Objective of Project 

Until now the mass transfer process is not fully understood. Observations also 
show that the orbital period of post-mass transfer systems are different and 
cannot be predicted by present models (see Figure [3}. Does there exist a corre- 
lation between those two aspects? The mass transfer induces a change in the 
angular momentum of both stars and might affect also the orbital period of the 
system. It will be a first step of the process of understanding mass transfer in 
binary systems by verifying the correlation between the barium abundance of 
a Barium star, being a binary star, and the orbital period of the system. Since 
barium is the most representative element of the mass transfer due to the en- 
hanced barium abundance, it seems reasonable to use it in order to verify the 
existence of a correlation. 



Figure 3: Graph showing the difference between the theoretical models and the 
observation of post-mass transfer stars 


(Izzard et al., 201071 


There were a few objectives to realize during this bachelor project. First, it 
was important to become familiar with binary systems and their interaction, as 
well as with stellar spectra. This latter theme has been understood by studying 
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the dependence of the stellar spectra on certain parameters. Once been in- 
troduced to these 'new' concepts, our investigation began with verifying, and 
even improving, the main stellar parameters of a sample of 9 Barium stars lo- 
cated at the main sequence of the H-R diagram accompanied by a white dwarf. 
Thereafter, the barium abundances were determined by investigating a few Ba 
II lines, that were thought to be representative for the determination of the bar¬ 
ium abundance. From these estimated abundances we sought for a correlation 
between them and their orbital periods, which were provided by our supervi¬ 
sor A. Escorza, who determined these using the radial velocity derived from 
their Doppler shift. 

It must be emphasized that during our research project a lot of assumptions 
were made, due to short period in which this project was realized. Ho wever, 
that does not imply that this research become less valuable; rather the objective 
was to obtain a qualitative idea on the correlation between the barium abun¬ 
dance and the orbital period than a precise and quantitative determination of 
the barium abundance. This correlation would have far-reaching implications 
in the understanding of the mass transfer process of binary systems. Therefore 
the importance of our research hides in the correlation and not in the precision 
of the estimated parameters of the Barium stars. 


2 Methods and Material 


2.1 HERMES and Atmospheric Models 

The observational data of our research project was provided by HERMES. This 
is one of the two main instruments of the Mercator telescope, which is a mod¬ 
ern semi-robotic telescope, installed at an observatory on La Palma, operated 
by the KU Leuven institute of Astronomy. The Mercator consists of two main 
instruments: HERMES and MAIA. Our main focus will be on HERMES (High 
Efficiency and Resolution Mercator Echelle Spectrograph), since it has two ad- 
vantages for this project. First, it allows us to study high-resolution astronomi- 
cal spectra, while on the other hand HERMES is already in operation for a few 
years which allows us to do measurements over relative long periods, which is 
the case for the Barium stars that are investigated (Raskin et al., 2011)| 


During the whole research the stellar spectral synthesis software PFANT 
written in Fortran language, was used to synthesize spectra. The MARCS- 
models of the mildly CN-cycle were used (|MARCS). This means that the soft¬ 
ware assumes the star of which the spectrum is synthesized to be spherical, 
in hydrostatic equilibrium and in LTE( Local Thermal Equilibrium). The LTE 
assumption simplifies a lot; the inflow of radiation is that of a black body with 
the same local temperature and to calculate the number of excited or ionized 
atoms the Boltzmann ([TJ and the Saha equations can be applied, resp. given 
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by 



(Eu ~ E x 

k B T 


( 1 ) 


with Ni the number of atoms in the lower energy level, N u the number of atoms 
in the upper energy level, gi and g u the statistical weight of both levels and Ei 
and E u the energies of the levels. 
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( 2 ) 


with N e the electron number density, m e the electron mass, N\ M the population 
in the ground state of two successive ionization stages of the atom and the 
ionization energy of level 1 (to reach the upper level u). 

The MARCS-modelling has also implications on the values of the different 
main parameters, T e ff ranges from 3500 K to 7000 K, logg from 0 to 3.5, and 
[M/ H] from —2 to 0.25. This means that PFANT can only synthesize the spec¬ 
tra with parameters within this ranges. 

2.2 Interpreting Astronomical Spectra 

Stars are characterized by several main parameters such as the luminosity (L), 
mass (M), radius (R), Chemical composition, etc. However, many of these main 
parameters are not directly measurable, and therefore other ways must be 
found to infer them. The two key parameters that can be observed are the 
effective temperature T e ff and the surface gravity, expressed by logg. These 
two observable parameters allow us to estimate the mass, luminosity and the 
radius of the star through the relations: 


L = inR 2 tjT* (( , 
GM 


(3) 


(4) 



where a is the Stefan-Boltzmann constant and G the gravitational constant. On 
the other hand, a secondary main parameter of a star is the stellar age, which 
can be estimated through evolutionary tracks in the H-R diagram |(Niemczura, 
jSmalley & Pych, 2014)| 

A few remarks on these relations and the observational parameters. First, 
it must be emphasized that these observational parameters only tells us some- 
thing about the surface of the star; the effective temperature is the temperature 
of a star if we assume the star to be a black body, while the surface gravity is the 
gravitational acceleration experienced at the surface of a star. A second thing to 
notice is that these relations reveal the interconnections between the effective 
temperature and the surface gravity. Hence, it must be emphasized that when 
estimating these two, it is highly recommended to solve them simultaneously. 
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Variation of the Effective Temperature 

The Balmer lines provide an excellent indicator to estimate the effective tem¬ 
perature for stars cooler than ±8000K (which is the case for the main sequence 
Barium stars studied) because they practically do not depend on the surface 
gravity However, some of these Balmer series are subjected to non-LTE effects 
because they are formed high in the stellar atmosphere, for example the core 
region of is strongly affected by NLTE effects. For that reason only two 
Balmer lines were studied to estimate the effective temperature, nl. and 
H 7 . To study these effects synthetic spectra of the Sun were created (see Figure 
|3]and|5]>. 

It can be seen from both hydrogen lines that as the effective temperature in- 
creases the wings of the spectral lines of and H 7 are broadened. This can de 
understood by the thermal Doppler effect. The motion of individual radiating 
atoms in a hot gas, responsible for line broadening, depend on the temperature 
of the gas in which they are moving. It can be found, using MaxwelTs velocity 
distribution, that this Doppler shift of the atoms (resulting in the broadening 
of the spectral lines) is given by 



(5) 


where k is the Boltzmann constant, and Ao the wavelength of the line centre 
|(Emerson, 1997)| 



Effect of Effective Temperature on H-beta line 
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Figure 4: The dependence of T e ff on the line of the Sun 
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Effect of Effective Temperature on H-gamma line 



Figure 5: The dependence of T e ff on the H 7 line of the Sun 


Since the Doppler shift is proportional to the width of the spectral lines, it 
follows that the width of the spectral lines is proportional to the effective tem¬ 
perature, i.e. if the effective temperature increases, the spectral lines become 
broadened, which is confirmed by the synthetic spectra created for the Sun. 

Variation of the Surface Gravity 

Although it has been remarked earlier that the effective temperature and the 
surface gravity are strongly related to each other, it is interesting to see how 
the variation of the surface gravity will affect the stellar spectra. To under- 
stand this dependence of the stellar spectra on the surface gravity, it is crucial 
to know what mechanisms are responsible for the broadening of the spectral 
lines. In fact, line broadening can de divided into 3 different categories: natural 
line broadening, which is due to Heisenberg's uncertainty principle, Doppler 
line broadening, due to the motion of the atoms in different directions and ve- 
locities, and pressure line broadening, due to effects of other particles on the 
radiating atom. The latter type of broadening will be of interest when studying 
the surface gravity. There are two important contributions to the pressure, nl. 
electron and gas pressure, denoted resp. by P e and P g . In Local Thermody- 
namic Equilibrium (LTE), the fraction of ions that can be absorbed is related to 
the electron pressure P e by the Saha equation j2|. 

From this equation it is clear that the line depth of certain ionized atoms 
depend on the electron pressure P e , which depends on the surface gravity. 
Hence, the surface gravity can be estimated by investigating weak lines that 
are formed by an element where most of it is in the next lower ionization stage. 
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for example Fe II will be a good indicator for the surface gravity, while Fe I will 
be insensitive to a change in surface gravity (see Figure[ó]and[7]). 

The inverse proportionality between the electron pressure and the line depth 
derived from the Saha equation (|2| can be confirmed by these two figures; as 
the surface gravity increases the line depth will decrease. Therefore, it can be 
concluded that these metal regions, where we have some good indicators such 
as Fe II will be a good measure for the surface gravity 


Variation of Surface Gravity 



Figure 6: The dependence of logg on Fe II lines of the Sun 
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Variation of Surface Gravity 



Wavelength (A) 

Figure 7: The dependence of logg on Fe II lines of the Sun 


Variation of the Metallicity 

A final main stellar parameter that may affect the stellar spectra is the metal¬ 
licity. The relation between the stellar spectra and metallicity is quite natural, 
since an increase in metallicity will result in an increase of metal lines and also 
an increase in their line depth. Consequently the metal regions used as a mea- 
sure for the surface gravity can also function as a measure for the metallicity 
of each star. Again to conclude this paragraph, it must be emphasized that all 
these three stellar parameters should be solved simultaneously. 


2.3 Determination Main Stellar Parameters 


Our first task of the project was to verify if the main stellar parameters from 
the literature (North et al., 1994 Takeda et al., 2005 Luck & Bond, 1991)] cor- 
respond to the observational data of HERMES. This was done by synthesizing 
a stellar spectrum based on the literature values of the star and to see if it ap- 
proximates the observational spectrum given by HERMES. For this purpose 
several regions of the stellar spectra has been investigated, in particular the 
and H 7 lines and the metal regions of the spectra (see discussion above). In the 
previous section it became already clear that these two hydrogen lines are quite 
sensitive to the effective temperature T e ff, while the spectral lines correspond- 
ing to the metals are sensitive to the surface gravity logg and the metallicity 
[M/H]. This method is illustrated in the figures |^,([9},([K)} and (??) below, 
where we used Barium star HD 89948 for which the literature values [(Luck &| 
|Bond / 1991) produce a synthetic spectrum that approximates the observational 
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data well. 


H-beta line of HD89948 



Figure 8: line of Barium star HD 89948 
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H-gamma line of HD89948 



Figure 9: H T line of Barium star HD 89948 


Furthermore, the literature values of the surface gravity and metallicity of 
the Barium star HD 89948 can be verified by studying the metal regions of the 
spectra. Again, it is found from figure [TÜJthat the synthetic spectrum produced 
by the literature values fit the observational HERMES data well, and thus it 
can be concluded that all stellar main parameters are accurate enough to allow 
further studies on the barium abundance of the star. 
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Metal region of Barium star HD89948 
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Figure 10: Metal region (5360,5390) A of Barium star HD 89948 


These first two figures ([8] and [9| show the and H 7 lines of the Barium 
star HD 89948 and it is clear that the synthetic spectrum fits the observational 
HERMES data quite well (neglecting the core), and thus it can be concluded 
that the effective temperature of the literature is a good estimate. 


Nonetheless, there were a few Barium stars for which the literature values 
did not produce a synthetic spectra that approximates the observational spec¬ 
tra well. Therefore an improvement of the main parameters was required for 
further study of these stars, and their barium abundance. This improvement 
of parameters was completed using the least-square method and a grid of pa¬ 
rameters. Since each parameter cannot be determined properly if the other 
parameters are kept constant, a grid was used to determine the main stellar 
parameters. To understand the concept, imagine a cube divided into smaller 
cubes (see Figure [TT]), where one axis represents the effective temperature, and 
the other two axes represent the surface gravity and metallicity, such that each 

point in the grid represents a set of parameters { T e ff, log g, 


}■ 


13 

































Figure 11: Visualisation of the concept of a grid 
IGridl 


This constmction permits to keep a continuous variation of parameters. To 
determine which set of parameters is the best approximation of the observa- 
tional parameters of the star, one can use the least-square method. This method 
computes the sum of all the squared differences of the synthetic and the ob- 
served spectrum for each wavelength, i.e. 


S = 


synthetic _ ypbserved 


( 6 ) 


where i runs over all wavelengths. The set of parameters where the least- 
square method leads to the smallest value S represents the best approximation 
of the observational main stellar parameters. 


2.4 Determination Barium Abundance 

For the barium abundance the same method can be used as for the determina¬ 
tion of the main parameters of a star, nl. the least-square method. To obtain a 
rough estimate of the abundance, a multiplot was generated with the GUI of 
PFANT for the values [Ba] = {—1,0,1,2,3} This was done for 8 barium lines 
(Ba II), for which 4 of them are resonant (see Table[I|. Then for each Ba II line 
a plot was made to verify if the Ba II line behaved well, or not, when other 
effects influenced them. In that case they are not representative/reliable to de¬ 
termine the barium abundance and are therefore neglected when determining 
the abundance. 

To concretize this selection of barium lines, it may be interesting to visual- 
ize a star for which a barium line was not representative for determining the 
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Table 1: The investigated Ba II lines 


Wavelength [A] 

Remarks 

4166.000 


4524.925 


4554.031 

Resonant 

4934.074 

Resonant 

5853.690 

Resonant 

6141.715 

Resonant 

6141.730 


6496.910 



abundance. This case is illustrated by the barium line at 5853.69A of Barium 
star BD+18 5215 (see Figure |Ï2|}. Although it might suggested that a higher 
barium abundance would fit better, this is not the case, since all for the other 
barium lines the barium abundances did produce a good fit for the HERMES- 
data. Therefore other effects are responsible for the broadening of this barium 
line. For example, it could be that there is a dependence on other parameters, 
or that NLTE effects are present. 


Barium line at 5853.69A 



Wavelength (A) 


Figure 12: barium line at 5853.69A of Barium star BD+18 5215 

After this selection of reliable Ba II lines, the best fit was determined for 
each line using the least-square method. Consequently, a second multiplot was 
made around the barium abundance corresponding to the best fit. In this way, a 
more precise estimate of the order 2 • 10 -2 was made for the barium abundance. 
Once an estimate was made for every reliable Ba II line, the mean of all the 
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abundances was calculated to find an estimate of the barium abundance of the 
star. This method was repeated for every star of the sample. 

2.5 Error Analysis 

Once the abundances and the orbital periods are determined, they can be plot- 
ted versus each other in a scatterplot together with their error-bars. This forms 
the last part of this section, namely determining the errors on the estimated 
barium abundances. 

If the total uncertainties on each abundance had to be computed, the follow- 
ing factors should be taken into account: the uncertainties on the parameters 
of the star, the relation between those parameters and the barium abundance 
together with the uncertainties due to the assumptions made (spherical model, 
LTE,...). Unfortunately, it is not realizable to compute this complete uncertainty 
of the abundances for each star. To obtain a global value of the abundance for 
each star, the mean of the established abundances for the different barium lines 
have been computed. However, computing a mean of a sample results in some 
statistical uncertainty over the mean of the total population, such that the mean 
of the population can be approximated by 

, s 

x = m ± —== 

Vn 

where x is the mean of the population (all the possible barium lines of the star), 
m the mean of the sample of barium lines for the star, s the Standard deviation 
of the sample, N the number of lines used in the sample. 

Hence, the statistical uncertainty is given by where s can be computed 
using 



For most of the stars the values for the main parameters were found in 
the literature. If in those papers the uncertainties of those stars were given, 
they also were adopted. When no uncertainties could be found, the gener al 
uncertainties (AT = ±10 0K, Alogg = ±1.0 dex and A [M/H\ = ±0.2 dex) 
were adopted. For the values established using the grid (star HD 76225) half 
of the length between two values of the grid were adopted as uncertainties on 
the multiple parameters. 


3 Results 

In this section the results of the different parts of the project will be presented. 
First the main parameters of the stars will be introduced, then the barium abun¬ 
dances and finally the plot of the orbital periods versus barium abundances. 
The sample used is: HD 22589, HD 15306, BD+18 5215, HD 34654, HD 95241, 
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HD 76225, HD 107574, HD 89948 and HD 150862. If the chosen stars were in 
the AGB phase, there would be other effects on the period. To avoid this, all 
the chosen Barium stars are on the main sequence in the H-R diagram. 

3.1 Main Parameters 

The literature values for the main stellar parameters of the Barium stars studied 
(see Table[2} produced synthetic spectra that were an acceptable approximation 
of the observational HERMES data, except for the Barium stars HD 15306 and 
HD 76225 (see Appendix A). The main stellar parameters were verified as dis- 
cussed in the section 'Determination Main Stellar Parameters', a plot was made 
of the H^ and H T lines for each star, where the wings of these lines were indi¬ 
cators of the effective temperature. On the other hand, a plot was made for the 
metal regions of each star, where we used Fe II lines and other ionized atoms 
as indicators of the surface gravity, while the metal lines in general served as 
indicators of the metallicity (see Appendix A). 


Table 2: Main stellar parameters of the sample Barium stars retrieved from the 
literature 


±100K, Alogg = ±1.0 dex and A[M/H] = ±0.2 dex 


Name 

T ef f [K] 

log g [dex] 

[M/H] [dex] 

HD 22589 

5500± 100 

2.30± 1.00 

-O.lOzb 0.20 

HD 15306 

6750± 100 

4.05± 1.00 

-0.34± 0.20 

BD+18 5215 

6290± 100 

4.49± 1.00 

-0.32± 0.20 

HD 34654 

6200 ± 100 

4.40 ± 1.00 

-O.lOzb 0.20 

HD 95241 

5904 ± 100 

3.70± 1.00 

-0.27zb 0.20 

HD 76225 

6010± 100 

3.97± 1.00 

-0.51zb 0.20 

HD 107574 

6250± 100 

2.90± 1.00 

-0.65zb 0.20 

HD 89948 

6000 ± 400 

4.00 ± 0.3 

-0.13 ± 0.05 

HD 150862 

6135± 100 

4.36± 1.00 

-0.18± 0.20 


However, two Barium stars (HD 76225 and HD 15306) of the sample did re- 
quire an improvement of the main stellar parameters. Despite the short period 
of time for this project, only for HD 76225 an improvement was made. From 
Figure[Ï3]and[T4]it is clear that the effective temperature from the literature is 
not a good approximation. Similarly, Figure [Ï5] also shows that an improve¬ 
ment of the surface gravity and/or metallicity is required. 
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H-beta line of HD76225 



Figure 13: line of Barium star HD 76225 


H-gamma line of HD76225 



Figure 14: H 7 line of Barium star HD 76225 
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Metal region of Barium star HD76225 



Figure 15: Metal region (5360,5390) A of Barium star HD 76225 

Therefore, a grid of five values for the effective temperature, five for the sur- 
face gravity and six for the metallicity was initially generated. Unfortunately 
using the MARCS-model of the mildly CN-cycle causes a small range of val¬ 
ues accessible for the parameters and most of the values for the surface gravity 
were outside this range. But PFANT offers the possibility to extrapolate points 
further than this range. However, af ter a quick examination of the extrapolated 
spectra, it was clear that the four values that feil out of the grid for the surface 
gravity were extrapolated to the same point. This meant that only two values 
for the surface gravity were available for the grid. 

All the spectra of the grid were analyzed over 5 different ranges: from the 

line to the telluric lines (4200 - 6700 A), around the line, around the H 7 
line, from the H 7 line to the the line and in a metallic region (4500 - 5400 
A). These 5 regions have been chosen to cut off most of the Non-Local Thermal 
Equilibrium effects on certain lines. It is also always better to investigate all the 
spectra over as much different regions as possible to get a realistic approxima- 
tion but it is recommended to avoid noisy regions and telluric lines (lines due 
to the Earth's atmosphere). For all of those regions the LS-method pointed dif¬ 
ferent spectra out of the sample but they all had the same temperature. It was 
then decided to work with this temperature and the metallicity found in the 
metallic region, because the other regions had values of the metallicity close 
to that of the metallic region. Then the LS-method ran over all spectra with 
constant effective temperature and constant metallicity and it pointed out one 
spectra that determined the value for the surface gravity. Hence, the best pa¬ 
rameters for HD 76225 are an effective temperature of 6200+2oo K, a log g of 
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3.97^4*4 dex (same as literature value), and a metallicity of (—0.30 ± 0.05) dex. 


- 0.34 

It was clear from Figure 


13 


that the line produced by the literature values 
did not approximate the observational HERMES data well. Now, it can be 
verified that the improvement of the main stellar parameters indeed solves the 
problem (see Figure |Ï6). 


H-beta line of HD76225 



Figure 16: line of Barium star HD 76225 after an improvement of the main 

stellar parameters 


3.2 Barium Abundances 

After all main stellar parameters had been verified or improved, the barium 
abundance could be determined following the procedure described in section 
'Determination Barium Abundance'. The selected barium lines that were con- 
sidered to be representative for the determination of the barium abundance can 
be found in Appendix B. 

From these selections of barium lines the best fit for each of them was deter¬ 
mined, using the least-square method. Consequently, the mean of all barium 
abundances corresponding to the best fit was calculated, resulting in the bar¬ 
ium abundances found in Table[3J where the errors on the barium abundances 
have been determined according to the method described in 'Error Analysis'. 
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Table 3: Barium abundances of the sample Barium stars 


Name Barium abundance [dex] 


HD 22589 
BD+18 5215 
HD 34654 
HD 95241 
HD 76225 
HD 107574 
HD 89948 
HD 150862 


(37 ± 9) • 10“ 2 
(87 ± 13) • ÏO" 2 
(17 ± 10) • 10“ 2 
(34 ± 8) • 10“ 2 
(128 ± 12) • 10“ 2 
(145 ± 6) • 10 -2 
(89 ± 9) • 10 -2 
(450 ± 96) • 10“ 3 


3.3 Barium Abundance versus Orbital Period 

Finally, the last step of this research project was to investigate if there is a cor- 
relation between the barium abundance of the Barium stars and their orbital 
period. The orbital periods (and the corresponding errors) of the sample of 
Barium stars, were provided by our supervisor A. Escorza (see Table |4]). 


Table 4: Orbital periods for the sample Barium stars 


Name 

Orbital period [days] 

HD 22589 

5518 ± 30 

BD+18 5215 

1517 ±15 

HD 34654 

976 ±1 

HD 95241 

5447 ± 30 

HD 76225 

2404 ± 7 

HD 107574 

1388 ± 13 

HD 89948 

667 ±2 

HD 150862 

290 ±2 


Combining the orbital periods (Table [2| and the determined barium abun¬ 
dances (Table [3} a scatter plot could be made (see Figure|T7|). From the scatter 
plot it is immediately clear that there is no linear correlation between both and 
thus no further statistical analysis have been proceeded. 
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Figure 17: Scatterplot of barium abundance versus orbital period 

4 Discussion 

As scientist it is important to be critical when analyzing and discussing your 
results, therefore a few remarks on the methods used and their reliability First 
of all, when investigating the main stellar parameters a quite mild verification 
of the literature values was performed. In the sense that no statistical methods 
were applied when verifying these literature values; when the synthetic spec¬ 
tra (produced by the literature values) fitted more or less the observational 
HERMES data, the naive assumption was made that these literatures values 
were adequate enough for further study However, a few of these literature 
values were outdated and thus it might be the case that they are not that ac¬ 
curate as was presumed. Some of them were also determined by photometry 
which is less accurate than the spectroscopie approach applied in this project. 
A better (also more time consuming) approach would be that a grid was gen- 
erated for each Barium star. To improve the established values by applying 
the LS-method on the grid, one should create another grid around the recently 
improved values and apply again the LS-method. By repeating this approach, 
the precision of the parameters will be increased. The values would also be 
improved if another software was used and other assumptions were made. 
Because of these assumptions, the PFANT-software has for example a certain 


Scatter Plot Banum abundance versus Orbital Period 
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range for the surface gravity, but unfortunately most of the literature values 
for the surface gravity were outside that range. Of course, there is still the pos- 
sibility to extrapolate to those values, but the uncertainty only increases and 
no good grid can be made around the values outside the range. It is always 
complicated to decide which approximations and assumptions must be made, 
because there are many factors that will affect a stellar spectrum and it is not 
always possible to know which ones are important. Unfortunately, all these 
uncertainties on the main stellar parameters persist when estimating the bar- 
ium abundances from the barium lines. 

Nevertheless, this does not imply that our results lose their value. In con- 
trary, their are a few important points that were taken into account that assure 
their reliability. For example, a careful treatment was made when studying 
the spectra, namely by taking NLTE effects into account. This was done when 
selecting the barium lines to estimate the barium abundance, as well as when 
investigating the main stellar parameters. Moreover, a restriction was made 
on the choice of Barium stars being studied; there were all chosen to be in the 
main-sequence, in order to avoid other external factors that may influence the 
observations. In addition, the data provided by HERMES certainly guarantees 
a high accuracy on the observational data. 

Since the stars from the sample are main-sequence stars, it should be ex- 
pected that no barium could be present and therefore no barium lines could be 
observed. Ho wever, from analysis it was verified that the stars from the inves- 
tigated sample do represent some barium lines. Therefore, it can be confirmed 
that the sample of stars being studied are indeed Barium stars. Nonetheless, 
one of the stars, nl. HD 34654 has an average barium abundance of (17 ± 
10) • 10 -2 dex, which is quite close to zero. This might indicate that HD 34654 
should not be part of the group of Barium stars, and thus further study is re- 
quired. The other abundances seem quite acceptable. There are also many 
other effects that have not be taken into account for the abundance of barium 
in a star. First, the abundances depend strongly on the stellar parameters. Sec- 
ondly, barium has different isotopes, and thus not the entire amount of barium 
in a star will be computed. Finally, there are the resonant lines where NLTE 
effects can occur. 

On the scatterplot it is clear that there is no correlation between the orbital 
period and the barium abundance, however it might as well be interesting to 
investigate other s-process elements (involving less quantum effects). This will 
provide a more definite answer on how the process of mass-transfer in binary 
Systems takes place. 


5 Conclusion 

During this project, we became familiar with studying stellar spectra; we rec- 
ognized their power of being able to study the physical aspects of objects that 
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are lightyears away from us, while on the other hand, we also realized their 
limitations. Our main objective was to obtain a better understanding of the 
mass-transfer process in binary systems by studying the correlation between 
the barium abundance of Barium stars and their orbital period. Unfortunately, 
thus far no conclusive resolution can be obtained about this interaction. It is 
clear that there is no (linear) correlation between the barium abundance and 
the orbital period. However, from a physical point of view one should expect 
that there is a correlation between the amount of mass transferred and their 
angular momenta (orbital period). Therefore we must infer that there are other 
parameters and external factors that play an important role in the process of 
mass-transfer. Secondly, this result indicates that barium is not an ideal ele¬ 
ment for studying this process, and one should rather (also) evaluate other s- 
process elements. On the other hand, it must emphasized that there also some 
limitations. For example, the companion of a Barium star, a white dwarf, does 
not allow a complete study of all its physical parameters. As a consequence, 
there are some hidden parameters in our study, for example we are not able to 
determine how much barium have been produced in the white dwarf. 

To summarize, despite the 'disappointing' result that there is no correlation, 
there are a still a lot of promising prospects for studying the process of mass- 
transfer in binary systems. To develop a better under standing, one should take 
into account more s-process elements, that are more reliable to study, while 
on the other hand a more careful determination of the main stellar parameters 
should be performed, in order to reduce the uncertainty that persist when es- 
timating the barium abundance from the barium lines. By taking into account 
these two improvements for studying the process of mass-transfer, one might 
take a first great step towards the understanding of this process and the inter¬ 
action of binary systems. 

To end this paper, we would like to thank our helpful supervisors Ana San¬ 
tos Escorza and Mike Laverick for all the advice and help during this bachelor 
project. Also a special thanks to Shreeya Shetye to help us providing good 
barium lines for studying the barium abundance, as well as to Dr. Andrew 
Tkachenko to warn us for NLTE effects when studying stellar spectra, and to 
suggest other regions for investigating the main stellar parameters. 
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Appendix A: The Main Stellar Parameters of Each 
Barium Star 

In this part of the appendix, the plots of the synthetic spectra approximating 
the observational HERMES data can be found. Moreover, these plots were 
made for the and H 7 line, and the metallicity regions of each barium star, 
to verify the main stellar parameters from the literature. 

HD 22589 

and H 7 lines 


H-beta line of Barium star HD22589 



Figure 18: line of Barium star HD 22589 


H-gamma line of Barium star HD22589 



Figure 19: H 7 line of Barium star HD 22589 
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Metallicity region 


Metal region of Barium star HD22589 



5375 5380 

Wavelength (A) 


Figure 20: Metal region (5360,5390) A of Barium star HD 22589 


HD 15306 

and H 7 lines 


H-beta line of HD15306 



Figure 21: line of Barium star HD 15306 
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H-gamma line of HD15306 



Figure 22: H 7 line of Barium star HD 15306 


Metallicity region 


Metal region of Barium star HD15306 



Figure 23: Metal region (5360,5390) A of Barium star HD 15306 


29 












BD+18 5215 

and H 7 lines 


H-beta line of BD18 5215 



Figure 24: line of Barium star BD+18 5215 


H-gamma line of BD18 5215 



Figure 25: H 7 line of Barium star BD+18 5215 
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Metallicity region 


Metal region of Barium star BD18 5215 



5370 5375 5380 

Wavelength (A) 


Figure 26: Metal region (5360,5390) A of Barium star BD+18 5215 


HD 34654 

and H 7 lines 


H-beta line of Barium star HD34654 



Figure 27: line of Barium star HD 34654 
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H-gamma line of Barium star HD34654 



Figure 28: H 7 line of Barium star HD 34654 


Metallicity region 


Metal region of Barium star HD34654 



5375 5380 

Wavelength (A) 


Figure 29: Metal region (5360,5390) A of Barium star HD 34654 
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HD 95241 


Hjg and H-y lines 


H-beta line of HD95241 



Figure 30: line of Barium star HD 95241 


H-gamma line of HD95241 



Figure 31: H T line of Barium star HD 95241 
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Metallicity region 


Metal region of Barium star HD95241 



Figure 32: Metal region (5360,5390) A of Barium star HD 95241 


HD 76225 

and H 7 lines 


H-beta line of HD76225 



Figure 33: line of Barium star HD 76225 


34 





















H-gamma line of HD76225 



Figure 34: H 7 line of Barium star HD 76225 


Metallicity region 


Metal region of Barium star HD76225 



Figure 35: Metal region (5360,5390) A of Barium star HD 76225 
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HD 76225 (Improved Main Stellar Parameters) 

Hjg and H-y lines 


H-beta line of HD76225 



Figure 36: line of Barium star HD 76225 


HD 107574 

and H 7 lines 


H-gamma line of Barium star HD107574 



Figure 37: line of Barium star HD 107574 
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H-gamma line of Barium star HD107574 



Figure 38: H 7 line of Barium star HD 107574 


Metallicity region 


Metal region of Barium star HD107574 



Figure 39: Metal region (5360,5390) A of Barium star HD 107574 
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HD 89948 


Hjg and H-y lines 


H-beta line of HD89948 



Figure 40: line of Barium star HD 89848 


H-gamma line of HD89948 



Figure 41: H T line of Barium star HD 89948 
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Metallicity region 


Metal region of Barium star HD89948 



Figure 42: Metal region (5360,5390) A of Barium star HD 89948 


HD 150862 

and H 7 lines 


H-beta line of HD150862 



Figure 43: line of Barium star HD 150862 
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H-gamma line of HD150862 



Figure 44: H 7 line of Barium star HD 150862 


Metallicity region 


Metal region of Barium star HD150862 



Figure 45: Metal region (5360,5390) A of Barium star HD 150862 
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Appendix B: The Ba II lines of Each Barium Star 

This part of the Appendix represents the multiplots of the barium lines that are 
considered to be reliable for the determination of the barium abundance. 

HD 22589 


Barium line at 41 66A 



Wavelength (A) 

Figure 46: Reliable barium line at 4166A of Barium star HD 22589 


Barium line at 4524.925A 



Figure 47: Reliable barium line at 4524.925A of Barium star HD 22589 
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Barium line at 4554.031A 



Figure 48: Reliable barium line at 4554.03A of Barium star HD 22589 


Barium line at 4934.074A 



Wavelength (A) 

Figure 49: Reliable barium line at 4934.07A of Barium star HD 22589 
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Barium line at 5853.69A 



Wavelength (A) 

Figure 50: Reliable barium line at 5853.68A of Barium star HD 22589 


Barium line at 6141.71 öA & 6141.73A 



Figure 51: Reliable barium line at 6141.715 A and 6141. 73A of Barium star HD 
22589 
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Barium line at 6496.91 A 



Figure 52: Reliable barium line at 6496.91A of Barium star HD 22589 
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BD+18 5215 


Barium line at 41 66A 



Wavelength (A) 

Figure 53: Reliable barium line at 4166A of Barium star BD+18 5215 


Barium line at 4554.031A 



Figure 54: Reliable barium line at 4554.03lA of Barium star BD+18 5215 
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Barium line at 4934.074A 



Figure 55: Reliable barium line at 4934.074A of Barium star BD+18 5215 

HD 34654 


Barium line at 41 66A 



Wavelength (A) 

Figure 56: Reliable barium line at 4166A of Barium star HD 34654 
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Barium line at 4524.925A 



Wavelength (A) 

Figure 57: Reliable barium line at 4524.925A of Barium star HD 34654 


Barium line at 4554.031A 



Figure 58: Reliable barium line at 4554.031A of Barium star HD 34654 
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Barium line at 4934.074A 



Figure 59: Reliable barium line at 4934.074A of Barium star HD 34654 


Barium line at 5853.69A 



Wavelength (A) 

Figure 60: Reliable barium line at 5853. 68A of Barium star HD 34654 
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Barium line at 6141.715A & 6141.73A 



Wavelength (A) 


Figure 61: Reliable barium line at 6141.715A and 6141.73A of Barium star HD 
34654 


Barium line at 6496.91 A 



Wavelength (A) 

Figure 62: Reliable barium line at 6469.91A of Barium star HD 34654 
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HD 95241 


Barium line at 4166A 



Figure 63: Reliable barium line at 4166A of Barium star HD 95241 


Barium line at 4524.925A 



Wavelength (A) 


Figure 64: Reliable barium line at 4524.925A of Barium star HD 95241 
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Barium line at 4554.031A 



Figure 65: Reliable barium line at 4554.031A of Barium star HD 95241 


Barium line at 4934.074A 



Figure 66: Reliable barium line at 4934.074A of Barium star HD 95241 
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Barium line at 5853.69A 



Wavelength (A) 


Figure 67: Reliable barium line at 5853. 6SA of Barium star HD 95241 


Barium line at 6141.715A & 6141.73A 



Wavelength (A) 


Figure 68: Reliable barium line at 6141.715A and 6141.73A of Barium star HD 
95241 
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Barium line at 6496.91 A 



Figure 69: Reliable barium line at 6469.91A of Barium star HD 95241 


HD 76225 


Barium line at 4554.031A 



Figure 70: Reliable barium line at 4554.031A of Barium star HD 76225 
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Barium line at 4934.074A 



Wavelength (A) 

Figure 71: Reliable barium line at 4934.074A of Barium star HD 76225 


Barium line at 5853.69A 



Figure 72: Reliable barium line at 5853. 68A of Barium star HD 76225 


54 















Barium line at 6141.715A & 6141.73A 



Figure 73: Reliable barium line at 6141.715A and 6141.73A of Barium star HD 
76225 


Barium line at 6496.91 A 



Wavelength (A) 

Figure 74: Reliable barium line at 6966.91A of Barium star HD 76225 
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HD 107574 


Barium line at 4524.925A 



Figure 75: Reliable barium line 4524.925A of Barium star HD 107574 


Barium line at 4554.031A 



Figure 76: Reliable barium line 4554.031 A of Barium star HD 107574 
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Barium line at 4934.074A 



Figure 77: Reliable barium line at 4934.074A of Barium star HD 107574 


Barium line at 6496.91 A 



Wavelength (A) 

Figure 78: Reliable barium line at 6496.91A of Barium star HD 107574 
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HD 89948 


Barium line at 41 66A 



Wavelength (A) 

Figure 79: Reliable barium line at 4166A of Barium star HD 89948 


Barium line at 4554.031A 



Figure 80: Reliable barium line at 4554.031A of Barium star HD 89948 
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Barium line at 4934.074A 



Wavelength (A) 

Figure 81: Reliable barium line at 4934.074A of Barium star HD 89948 


Barium line at 5853.69A 



Wavelength (A) 

Figure 82: Reliable barium line at 5853. 68A of Barium star HD 89948 
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Barium line at 6141.715A & 6141.73A 



Wavelength (A) 


Figure 83: Reliable barium line at 6141.715A and 6141.73A of Barium star HD 
89948 


Barium line at 6496.91 A 



Figure 84: Reliable barium line at 6496.91A of Barium star F1D 89948 
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HD 150862 


Barium line at 41 66A 



Wavelength (A) 

Figure 85: Reliable barium line at 4166A of Barium star HD 150862 


Barium line at 4524.925A 



Figure 86: Reliable barium line at 4524.925A of Barium star HD 150862 
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Barium line at 4554.031A 



Figure 87: Reliable barium line at 4554.031A of Barium star HD 150862 


Barium line at 4554.031A 



Wavelength (A) 

Figure 88: Reliable barium line at 4934.074A of Barium star HD 150862 
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Barium Abundances Estimated from Barium lines 


Table 5: Estimated barium abundances for the reliable barium lines 


Barium lines at 


4166A 4524.925A 4554.031A 4934.074A 5853.69A 6141.715A 6496.91A 


HD 22589 

0 

0.4 

0.2 

0.2 

0.6 

0.6 

0.6 

BD+18 5215 

0.6 

NA 

1 

1 

NA 

NA 

NA 

HD 34654 

-0.2 

0.2 

0.2 

0.4 

-0.2 

0.4 

0.4 

HD 95241 

0 

0.2 

0.2 

0.4 

0.4 

0.6 

0.6 

HD 76225 

NA 

NA 

1 

1 

1.6 

1.4 

1.4 

HD 107574 

NA 

1.4 

1.4 

1.4 

NA 

NA 

1.6 

HD 89948 

0.8 

0.6 

0.8 

0.8 

1 

1.2 

1 

HD 150862 

0.2 

0.4 

0.6 

0.6 

NA 

NA 

NA 
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